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Previous studies indicated that intrinsic point defects play a crucial role for the density of states 
of ferromagnetic half-metals in the band gap region: at large concentrations, defect-derived bands 
might close the gap at the Fermi energy in the minority spin channel. In this work, structural 
disorder in the Co- and Mn-sublattices of the full Heusler alloy Co2- :E Mni+ a; Si ( — 1 < x < 2) is 
investigated with a cluster expansion approach, parametrized using all-electron density functional 
theory calculations. By establishing two separate cluster expansions, one for the formation energy 
and one for the total spin moment, we are in a position to determine the stability of different 
configurations, to predict new (also half-metallic) ground states and to extend the known Slater- 
Pauling rule for ideally stoichiometric Heusler alloys to non-stoichiometric, Mn-rich compositions. 
This enables us to identify potentially half-metallic structures in the Mn-rich region. With the 
help of Monte Carlo simulations based on the cluster expansion, we establish theoretically that 
Co2- a; Mni+ :E Si close to the stoichiometric composition ought to show a high degree of structural 
order in thermodynamic equilibrium. Hence, samples prepared with the correct stoichiometry should 
indeed be half-metallic after thermal annealing. Moreover, we predict that adding a small amount of 
Mn to stoichiometric Co2MnSi allows suppression of the thermally activated formation of detrimental 
Co antisites. At Mn-rich compositions (x > 1), the ordered ground state structures predicted for 
zero temperature are found to be thermally unstable and to decompose into Co2MnSi and MnaSi 
above room temperature. 



PACS numbers: 71.20.Be, 75.50.Cc 

I. INTRODUCTION 

Since their prediction in 1983 (1) half-metallic ferro- 
magnets have attracted much interest from a fundamen- 
tal scientific point of view and as materials for spintron- 
ics devices. But there has been and still is a strug- 
gle to verify the existence of the gap in the minority 
spin channel experimentally. Following one line of re- 
search, the difficulty to observe a gap at finite temper- 
atures is an intrinsic property of a ferromagnetic half- 
metal: In a many-particle description, thermal excita- 
tions of the electronic system formed by superpositions 
of majority spin states and virtual magnons appear in the 
Kohn-Sham band gap. These states filling up the half- 
metallic gap are referred to as non-quasiparticle states^. 
Other researches have invoked extrinsic factors, such as 
structural defects, including atomic disorder, surface or 
interface states, to explain the low experimental spin 
polarization4££il For half-metallic Heusler alloys first- 
principles calculations for bulk systems have shown that 
disorder in the sublattices 8 - or impurity bands induced by 
a high concentration of point defects^ may close the gap. 

Here, we are interested in the full Heusler alloy 
Co2MnSi because it is a candidate for spintronics devices 
working at room temperature due to its large (predicted) 
Kohn-Sham gap ; 10 ' 11 and its high measured Curie tem- 
perature Tc = 985 K i 12 i 13 On the other hand, experimen- 
tal values for the spin polarization of only P = 50 — 60% 



(instead of the expected 100%) have been measured for 
crystalline bulk and film samples using point-contact An- 
dreev reflexionj^i 7 - Somewhat higher values in the range 
of P = 60 — 90% have been inferred from the tunnel- 
ing magneto resistance ratio measured in tunnel junc- 
tions with a Co 2 MnSi electrode&i^±&±£ at low tem- 
peratures, however, with a significant drop of P around 
room temperature . 18 ' 19 Already in the early works, in- 
sufficiently ordered samples, in particular due to Co-Mn 
swaps, have been made responsible for the low values 
of P. Neutron diffraction measurements show that 14% 
of the Mn sites are occupied with Co atoms in some 
samples.— Density functional theory calculations using 
supercell geometries for Co2MnSi and Co2MnGe show 
that some defects may change the electronic structure 
drastically; 9 ^' but only four types of defects have been 
investigated so far. 

The aim of this work is to study realistic, partially 
disordered Co2MnSi with ab initio accuracy over a wide 
concentration range and at elevated temperatures. We 
are especially interested in non-stoichiometric, Mn-rich 
Co2MnSi because the results of single defect calculations 
suggest that additional Mn atoms at Co sites may help 
to preserve the band gap.— Moreover, disorder, "frozen 
in" during sample preparation, could help to rationalize 
the non-ideal behavior of Co2MnSi observed in many ex- 
periments. In this study, we aim at identifying a range 
of compositions and annealing temperatures where half- 
metallicity can be expected. 
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This paper is organized as follows: First, we briefly 
describe the cluster expansion method based on density 
functional theory calculations. It enables us to explore 
formation energies and magnetism in a huge configura- 
tion space. In Section 3, the properties of isolated defects 
and of the ground-state structures at zero temperature 
are reported. Moreover, we propose an extension of the 
Slater-Pauling rule that is useful in assessing the mag- 
netic properties of Mn-rich Co2MnSi. Finally, the cluster 
expansion of the formation energy is employed to perform 
Monte Carlo simulations both for nearly ideal Co 2 MnSi 
at finite temperatures, and for a newly discovered or- 
dered alloy of composition Co2Mii4Si2. In Section 4, we 
summarize our results and conclude. 



II. METHODS 

A. The cluster expansion 

The cluster expansion (CE) method can be used to 
parametrize any physical property F of a multicom- 
ponent system that depends uniquely on the atomic 
configuration^ An excellent review of the method and 
its applications is given in Ref. [23|. Here, we only give a 
brief introduction of the underlying principle. 

Let us consider a binary alloy Ai-^B^ with N lattice 
sites. Each site is described by an occupation variable a 
that is +1(— 1) if a lattice point is occupied by atom A 
(B). A configuration of the whole crystal is then charac- 
terized by the occupation vector er = {<ti, 02, . . . , ctjv}- 
The property F can be expanded into a series with a 
structure similar to the Hamiltonian of the Ising model 

F((t) = J + ^ Ji<Ji +^ J, ,rr,rr , + ./, ,;. n,a ,n h + . . . 

i i.j 

where the pairs, triplets and higher order terms are the 
so-called "figures" of the CE and the Jj are the effective 
cluster interactions (ECI). Equation (JTJ) is in principle 
exact*22 but has to be truncated for practical reasons. 
The unknown parameters Jj of the CE are determined 
from the results of a relatively small number of configu- 
rations (1=3 50) obtained through first-principles computa- 
tions using a procedure that minimizes the mean-square 
deviations between the predicted and the calculated val- 
ues of F(tr) on subsets of the structures that have been 
calculated within DFTi^i This approach is also called the 
structure inversion method^ There remains the crucial 
question of how to identify the physically relevant in- 
teractions. Minimizing the least-squares error of the fit 
is not sufficient, since it does not guarantee the trans- 
ferability of the fit to other (unknown) configurations. 
A more suitable measure is the cross-validation score 
(CVS). Here, we use the leave-one-out cross-validation 



score, defined by 

1 n 

C 2 cvs = -Y,( F r~FC?)\ (2) 

i 

where the F^ E are predicted by a least-squares fit to n— 1 

first-principles (FP) values excluding Ff p . Therefore, 
Ccvs gives an estimate of the predictive power of the 
fit. The set of figures that minimizes the cross-validation 
score is termed the optimal cluster expansion. 

The method can be extended to multi-component 
systems.— The CE Hamiltonian is most widely used for 
the calculation of alloy formation energies and as input 
for Monte Carlo simulations of thermodynamic quanti- 
ties and phase diagrams, e.g. Cu-Au intermetallics^ or 
concentration of vacancies in the Al^^Li^. alloy. 27 The 
method has also successfully been applied for the calcu- 
lation of magnetic properties of Fe-Co alloys 28 - and the 
computation of Curie temperatures of Mn-doped GaAs. 29 

B. Computational details 

Density functional theory (DFT), using the highly 
accurate all-electron full-potential linearized augmented 
plane wave (FP-LAPW) method implemented in the 
Wien2k code^S, has been used to calculate total ener- 
gies and magnetic spin moments of different composi- 
tions of Coa;Mn a Si z . The generalized gradient approxi- 
mation (PBE962I) to the exchange-correlation functional 
has been adopted because the GGA gives better results 
than the LDA (regarding lattice constants and bulk mod- 
uli) for Co 2 MnSi according to Refs.fTll and[32l. For inves- 
tigations of atomic defects, 2 x 2 x 2 fee supercells with 
32 atoms and a 6 x 6 x 6 k-point mesh have been used. 
As input for the cluster expansion, unit cells containing 
4 to 32 atoms have been considered with dense k-point 
meshes (12 x 12 x 12 for 4 atom cells, 6 x 6 x 6 for 32 atom 
cells) to ensure convergence. The differences caused by 
the different k-point meshes in different unit cell shapes 
are very small (e.g. the total energy difference between 
Co2MnSi in a 4-atom fee unit cell and a 16-atom sim- 
ple cubic unit cell is less than 0.1 me V/ atom) and do 
not influence the results. The muffin-tin radii were set 
to R MT (C6) = R MT (Mn) = 2.1 bohr and R M t(Si) = 
2.0 bohr. The spin magnetic moments per atom quoted 
below have been calculated by integrating the spin den- 
sities obtained from converged self-consistent calcula- 
tions within the muffin-tin spheres, using these radii. 
Larger muffin-tin radii (touching spheres are achieved 
with R MT (Co) = R M t{Mii) = R M t($i) = 2.3 bohr) 
have been tested as well, but are unsuitable for the re- 
laxation of the structures while they have negligible ef- 
fect (less than 3%) on the analysis of the atom-resolved 
spin moments. In the interstitial region a plane wave ex- 
pansion cutoff of 15.5 Ry was used. All internal atomic 
coordinates were relaxed until the forces on the nuclei 
were smaller than 2 mRy/bohr. 
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FIG. 1: (color online) The L2i structure of Co2MnSi. Green 
(large) - Co atoms, pink (medium) - Mn atoms and blue 
(small) - Si atoms. 



III. RESULTS 

A. Ideal Co^MnSi and the role of defects 

The L2i structure (common for all full Heusler alloys) 
of Co2MnSi is shown in Fig. [TJ Within density func- 
tional theory (both LDA and GGA) the half-metallicity 
of ideal Co2MnSi is well established and a summary of 
some calculated properties is given in Table [T] Through- 
out the paper we use the term "spin gap" for the energy 
gap between the highest occupied state in the minority 
spin channel and the Fermi energy and the term "band 
gap" for the energy gap between the highest occupied 
state and the lowest unoccupied state in the minority 
spin channel. 

A scheme for the hybridization of the o?-bands in the 
full Heusler alloys that leads to the formation of the mi- 
nority band gap has been established by Galanakis et al. 
(Ref. [33h . Most of the Heuslers follow a Slater-Pauling 
(SP) rule that relates the total spin moment M to t of a 
compound to its number of valence electrons Ny as fol- 
lows: 

M tot = N V -2A. (3) 

In other words, the SP rule states that the number of elec- 
trons in the minority spin channel is always 12 since this 
corresponds to a filling of four sp bands, two ti g bands 
(occupied by three electrons each) and one e g band (oc- 
cupied by two electrons) in the spin-down channel. These 
filled bands are separated from the empty bands by a hy- 
bridization gap. The difference in the number of valence 
electrons due to introducing other elements is accommo- 
dated in the majority spin channel only. A necessary (but 
not sufficient) condition for half-metallicity is an integer 
total spin moment per unit cell. The SP rule formulates 
a stricter, material specific (but still not sufficient) con- 
dition for the half-metallicity of Heusler alloys. 

Next we report results for several kinds of defects. The 
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a (A) 


5.629 


5.634 


exp. 


5.607 


E s (eV) 


0.35 


0.33 


n/a 


0.36 


En feVl 


0.84 


0.81 


n/a 


0.42 


mco (/xs) 


1.07 


1.06 


1.02 


n/a 


m Mn (mb) 


2.91 


2.92 


2.97 


n/a 


m S i (a*s) 


-0.04 


-0.04 


-0.07 


n/a 


Mtot (/is) 


5.00 


5.00 


5.00 


n/a 



TABLE I: Calculated (DFT) lattice constant a, Kohn-Sham 
spin gap Es, Kohn-Sham band gap Eg, atomic spin moments 
rrii (i — Co, Mn and Si) and total spin moment per unit cell 
Mtot of bulk Co 2 MnSi. 

defect formation energy is calculated as 34 

E f (X) = E tot {X) - £ tot (Co 2 MnSi) - ^n ilH (4) 

i 

where E to t(X) is the DFT total energy of the supercell 
containing defect X and £t ot (Co2MnSi) is the DFT to- 
tal energy of ideal Co2MnSi in an equivalent supercell. 
The difference in the number of atoms to the stoichio- 
metric composition is taken into account by rii that is 
+1(— 1) for an excess (deficiency) of atom of species i, 
and is the chemical potential of the reservoir (total 
energy per unit cell). We calculate the formation ener- 
gies with respect to two different sets of reservoirs. The 
set A of reservoirs is the "conventional" one where we use 
the elements Co, Mn or Si in their most stable bulk form, 
i.e. ferromagnetic hep Co, cv-Mn (calculated as antiferro- 
magnetic 7-Mn plus a correction for the energy difference 
between a- and 7-Mn taken from Ref. 1351 ) and diamond 
Si. Since we are especially interested in the Co-Mn in- 
teractions in a Co- and Mn-rich regime, no elemental Si 
will be present, but compounds such as Mn 3 Si or CoaSi 
may be formed instead. Therefore, we define a new set 
of reservoirs B where we eliminate fisi from Eq. ([¥]) us- 
ing the laws of mass action for formation of M^Si and 
C03SL What enters in Eq. (0J are now the total ener- 
gies of hep Co, a-Mn, CosSi and M^Si (all in their bulk 
form). The formation energies with respect to both sets 
of reservoirs, magnetic moments and spin gaps are shown 
in Table [Til For the cases that have already been studied 
in Ref. |9| we find excellent agreement for the magnetic 
and electronic properties. 

For most defective supercells the magnetic moments 
of atoms at regular sites are barely changed. The two 
exceptions are Si atoms at Co sites and the vacancies. Si 
atoms at Co sites increase the spin moments of the Co 
atoms surrounding them to 1.2 /ig. In the case of the Sico 
antisite this leads to a total spin moment of 41 instead 
of the expected 39 hb- In case of a vacancy, the moments 
of the neighboring atoms are reduced, both for Co (0.8 — 
0.9 hb) and Mn (2.7 — 2.8 /is) neighbors. A similar 
behavior is observed for Si, Mn and Co vacancies. The 
defects themselves show different behavior. Co atoms at 
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Mn sites have magnetic moments of 1 fiB (Co-Mn swap) 
and 1.2 /is (CoMn antisite). Co atoms at Si sites show 
an increased magnetic moment of about 1.6 fiB (for swap 
and antisite). In agreement with earlier calculations 36 , 
we find that Mn atoms at Co sites (as part of a swap and 
as antisite) couple antiparallel to atoms at regular sites 
with a magnetic moment of around — 1 fiB which leads to 
reduced total magnetic moments compared to the ideal 
case. 

For a Mnco antisite in a 32-atom supercell we have 
done a calculation with a fixed total spin moment (FSM) 
of the supercell of 42 \ib which corresponds to a substi- 
tution of niQ s» 1 Hb with tomh ~ 3 fig. In this case 
the Mn antisite atom changes its magnetic moment from 
— 1 fig to 2.3 fiB and the Co atoms (at regular sites) 
show slightly increased spin moments of about 1.15 fiB- 
The energy difference between the FSM calculation and 
the free one with Mtot = 38 fiB is about 1.1 eV. This 
reflects that the antiparallel coupling of the Mn atom at 
a Co site is strongly preferred over parallel coupling. 

As an example for the impact on the electronic struc- 
ture the total density of states (DOS) of the Co-Si swap 
is compared to the ideal case in Fig. [2] There are two 
defect states close to the Fermi energy that are solely 
derived from the Co atom at Si site. One of them is 
occupied and lies slightly below Ep, the other one is un- 
occupied. Due to these states the remaining spin gap 
is very small. The DOS of the Cosi antisite is similar 
to that of the Co-Si swap showing also the two defect 
states. But in this case the occupied state lies at the 
Fermi energy and closes the spin gap. 

We note that the described position of levels with re- 
spect to the Fermi energy refers to the Kohn-Sham en- 
ergy levels calculated within the DFT-GGA approxima- 
tion. We expect the qualitative picture, e.g. the orbital 
character of the defect states, to remain correct if meth- 
ods beyond DFT-GGA were used. The relative energy 
position of Mn and Co orbitals is correctly described by 
DFT-GGA, and hence the predicted Co orbital charac- 
ter of highest occupied spin-down states will persist even 
in a higher-level method. However, the exact position of 
the Fermi level relative to the defect-derived states may 
depend on the treatment of electronic exchange and cor- 
relation effects. For example, in recent calculations per- 
formed with the GGA+U method it was observed that 
the position of the Fermi level can be shifted relative to 
the band edges in the spin down channel if the value of 
U is varied^ 

Taking into account our results together with the 
known results for Co-Mn swaps, CoMn, Mnc antisites^, 
Mnsi and SiMn antisites^ and vacancies 3 ^ (all obtained 
in the framework of DFT) we can conclude that defects 
with Co atoms at Mn or Si sites may close the spin gap 
while all other defects have only minor effects on the elec- 
tronic structure and are harmless for the half-metallicity. 
However, the vacancies lead to a considerable reduction 
in the spin gap (see Table HI)) . and have been predicted 39 
to disturb the half-metallicity for concentrations > 10%. 





Ef (eV) 


Ef (eV) 
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E s (eV) 


CoMn antisite 


0.75 


0.65 


38.0 


0.05 


Cosi antisite 


2.25 


1.26 


41.0 


0.00 


Mnco antisite 


0.36 


0.45 


38.0 


0.39 


Mnsi antisite 


1.69 


0.40 


43.0 


0.42 


Sico antisite 


1.78 


2.77 


41.0 


0.59 


SiMn antisite 


-0.21 


1.08 


37.0 


0.33 


Co-Mn swap 


1.20 


1.20 


36.0 


0.20 


Co-Si swap 


3.84 


3.84 


42.0 


0.12 


Mn-Si swap 


1.17 


1.17 


40.0 


0.43 


Vco vacancy 


1.04 


1.04 


37.0 


0.26 


Vmu vacancy 


1.43 


1.43 


33.0 


0.15 


Vsi vacancy 


3.74 


2.72 


36.0 


0.14 



TABLE II: Formation energies Ef according to eq. @ with 
respect to reservoirs A (Co, Mn and Si) and B (Co, Mn, CoaSi 
and MnaSi), total spin moments Mtot and Kohn-Sham spin 
gap Es of some possible defects in supercells with 32 atoms. 
Stoichiometric Co2MnSi in a 32-atom supercell has a total 
spin moment of 40 fiB- 
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FIG. 2: (color online) The total density of states of ideal 
Co2MnSi (black dotted line) and the Co-Si swap (red line). 
The crucial difference lies in the states in the band gap , i.e. 
between — 0.5eV and +0.35eV, indicated by the two arrows, 
that are derived solely from the Co atoms at Si sites. 



B. CE for pseudo-binary C02 ^Mni+^Si 

The results for single atomic defects show that the elec- 
tronic and magnetic properties of Co2MnSi are strongly 
altered by defects involving Co (leading to defect states 
close to or even at the Fermi energy) and Mn (magnetic 
moment antiparallel to regular Co and Mn) atoms. We 
verified by explicit calculations that other defects, such 
as vacancies, Mn-Si swaps or Si antisites have less impact 
on the electronic structure (i.e. they do not induce de- 
fect states in the half-metallic gap) and on the magnetic 
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FIG. 3: (color online) Optimal CE figures labeled by their in- 
dex. The green (large) circles are Co sites, the pink (medium) 
circles Mn sites, and the blue (small) circles Si sites. The 
triples 14 and 15 (not shown) are the smallest possible equi- 
lateral triangles in Co(lll) and Mn(lll) planes. Co-Mn pair 
5 and Co-Co-Mn triple 10 lie outside the lattice planes shown 
in the figure. 



properties (no antiparallel coupling) . Since the formation 
energies of a Co-Si swap and a Cosi antisite are much 
larger than the corresponding Co-Mn swap and CoMn 
antisite (for both sets of reservoirs considered in Tab. 
HI)) , the former two are neglected in the following. Thus, 
we regard the material as a pseudo-binary alloy where 
only the interactions between the Co and Mn atoms are 
taken into account while the Si sublattice is kept fixed. 

The cluster expansions are constructed using the 
ATAT software3il This 'alloy theoretic automated 
toolkit' generates structures on the given parent lattice 
and uses the DFT values (e.g. total energy) of these 
structures to fit all interactions up to a given maximum 
number of points of the CE figures (here six). It com- 
pares the cross-validation scores of the different CEs and 
chooses the CE with the lowest one. 

In the following the CE for the formation energies of 
the compounds is called ECE and the one for the total 
spin moment is called MCE. The total energies and spin 
moments of 60 structures have been computed with DFT- 
GGA for the fit of the interactions. The optimal set of fig- 
ures for both CEs consists of 9 pairs and 6 triplets (most 
of them shown in Fig. [3| with cross validation scores 
C^f = 5.8 meV/atom and = 0.05 ^ s /atom. As 

will be seen later, Cqy§ is of the same order of magni- 
tude as the ordering energies of Mn-rich structures. The 
strengths of the effective cluster interactions (ECI) are 
shown in Fig. 0] both for the CE of the total energy and 
the CE of the total spin moment. In the energy expan- 
sion, the leading attractive interactions (negative values) 
are the pairs 2 and 4 and the triple 14. These interactions 
ensure that the sites in the Co(100) plane, the Mn sites in 
the MnSi(100) plane, and the sites in the Co(lll) plane, 
respectively, are preferentially occupied by like species of 
the suitable kind. Tests with a CE expansion containing 
more long-ranged pair interactions, or additional triples 
or quadruples, didn't lead to a significant reduction of 

riECE 

^cvs ■ 

One has to keep in mind that Co 3 Si and Mn 3 Si are the 
reference structures and reservoirs for the ECE and that 
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FIG. 4: Effective cluster interactions (ECI) for the CE for 
the formation energy (top), and for the total spin moment 
(bottom). For the indexing of the CE figures, cf. Fig. [3] 



the Si sublattice is fixed. Therefore, the general eq. ((3]) 
is adopted for this special case, by eliminating //gi in the 
way described above. Every structure contains exactly 
25% Si and the chemical potentials are now the total 
energies per unit cell of Co 3 Si and Mn 3 Si. The formation 
energy per atom of a Co2- :r Mni_|_ ;r Si compound is then 
defined as 



fiECE 

1 



f (Co a Mn b Si( a+fc) / 3 ) = 

E(a,b) - |s(Co 3 Si) - ^(Mn 3 Si) 



a 



(5) 



where E(a, b) is the total energy of the structure calcu- 
lated with DFT. 

In Fig. [5J the values from the DFT-GGA calculations 
are compared to the CE values. Excellent agreement is 
found. Additionally, a ground state search has been per- 
formed and the ground state line is included in the up- 
per part of Fig. [5j For three structures a, 7 and (3 with 
concentrations x(a) < x(j) < x{[3) that have the low- 
est formation energies at their respective concentrations 
structure 7 has to fulfill the condition 



E ECE (^\ < ^(7) X \P) gECE 

* x(a) — x((3) ' 



(a)- 



#(7) — x(a) 
x{j3) — x(a) 



E 



ECE 



( 6 ) 

to be part of the ground state line at £(7). 23 In other 
words, the ground state line is the lower convex enve- 
lope of all data points in the formation energy diagram. 
With the ECE Hamiltonian the formation energies of all 
combinatorial possibilities in a unit cell with 32 atoms 
(24 Co/Mn sites, 2 24 w 16.8 million) and of 10 7 struc- 
tures in a unit cell with 128 atoms (96 Co/Mn sites) 
where the sites have been randomly occupied with Co 
or Mn have been computed. The well-known L2i struc- 
ture of Co2MnSi is the most stable configuration. Four 
new ordered ground states at zero temperature are pre- 
dicted at Mn concentrations of x = 0.1875, x = 0.5, 
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FIG. 5: Comparison of the values from DFT-GGA and from 
the CEs (top: formation energy, bottom: total spin moment). 
The excellent agreement between the the values demonstrates 
the quality of the two fits. 

x = 0.625, and x = 0.6875. At these concentrations, 
slight (almost invisible) cusps appear, deviating by about 
5 mV/atom from the convex envelope. Since these de- 
viations are very small, test calculations have been per- 
formed, demonstrating that the formation energies pre- 
dicted by the ECE agree with the DFT-GGA values (cf. 
Fig. [5]). For concentrations different from the ordered 
ground states, the lowest energy state of the system at 
zero temperature consists of two regions of coexisting 
phases. The fact that the data points in Fig. [8] do not 
reach down to the convex envelope at all concentrations 
is due to the finite size of the unit cell used, which doesn't 
allow for phase separation to occur. 

C. The ordered ground states 

The ordered ground state at the Co-rich side is a sin- 
gle CoMn antisite in a 16-atom simple cubic unit cell 
(CogMn3Si4). ft has a broad Co-derived band crossing 
the Fermi energy, leading to a small spin polarization 
and rendering this structure uninteresting for potential 
spintronics applications. 

We will concentrate on the ordered ground states at 
the Mn-rich side. They have highly symmetric struc- 
tures that show a certain trend in the incorporation of 
the Mn atoms at Co sites as shown in Fig. [5] First, in 
Co2Mn4Si2 each second layer of Co atoms is completely 
replaced by Mn atoms. Increasing the Mn concentra- 



FIG. 6: (color online) Structures of the Mn-rich ground states, 
top left: Co2Mn4Si2, bottom left: Co2MnioSi4 and right: 
CoMnnSi4. Green (large) - Co atoms, pink (medium) - Mn 
atoms, and blue (small) - Si atoms. 

tion further leads to additional Mn rows in the remaining 
Co planes (Co2Mni Si4 and CoMnnSi4). These ground 
states follow the general trend in the physics of alloys 
that stable structures have relatively small and highly 
symmetric unit cells. 

Ideal Co2MnSi is a half metal according to DFT cal- 
culations. The electronic structures of the other ground 
states show one remarkable result. The total DOS of 
Co2Mn4Si2 has a gap in the minority spin channel around 
Ep and is therefore a half metal, too. Investigating the 
orbital-resolved DOS (shown in Fig. [7J reveals that the 
spin-down DOS of the Mn atom at Co site is almost iden- 
tical to the spin-down DOS of Co. This behavior is also 
reflected in the magnetic moments of the Mn atoms at 
Co sites. They align antiparallel to the regular Co and 
Mn magnetic moments, and their magnitude ranges from 
—0.6 fin to —0.7 (j.b- The integration of the DOS from 
-co to Ep confirms that the missing electrons (due to 
replacing Co with Mn) in this structure are taken only 
from the majority spin channel. This complies with the 
Slater-Pauling rule for the full Heusler alloys mentioned 
earlier. 



D. Extended Slater-Pauling rule 

In the supercell calculations, Mnc antisites reduce the 
spin moment by 2 \ib (from +1 fJ-B to —I /is) which 
corresponds exactly to the reduction of the number of 
valence electrons (from 7 Co d electrons to 5 Mn d elec- 
trons). This suggests that the Slater-Pauling rule may 
still be valid even for non-stoichiometric, Mn-rich compo- 
sitions Co2- x Mni +a .Si. We have performed DFT-GGA 
calculations of Mn-rich structures to verify the linear re- 
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FIG. 7: (color online) The density of states projected on the d- 
orbitals of the predicted half-metallic ground state Co2Mn4Si2 
(black line - Co, blue dash-dotted line - regular Mn, and red 
dashed line - Mn atom at Co site). In the minority spin 
channel the DOS of the Co atom and of the Mn atom at the 
Co site are almost identical. 



FIG. 8: (color online) Direct enumeration of the formation 
energies according to eq. ([5]) of 27 million structures. The 
black dots belong to structures that obey the extended Slater- 
Pauling rule of eq. For the structures marked by blue 
circles it has been confirmed by DFT-GGA calculations that 
they have a spin gap and and satisfy eq. 0. 
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96 
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83 
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72 


0.19 
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54 
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72 
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96 
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Co 2 Mni Si 4 
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8.0 


96 
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Co 2 Mni Si 4 
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0.04 


CoMn 8 Si 3 


77 


5.0 


72 


0.00 


CoMnnSu 


102 


6.02 


95.98 


0.00 


CoMnnSu 


102 


6.02 


95.98 


0.00 


CoMnnSu 


102 


6.01 


95.99 


0.00 



TABLE III: Number of valence electrons Nv, total spin mo- 
ment Mtot, and spin gap Es of some Mn-rich structures. 



The relation between the number of valence electrons 
and the total spin moment can be expressed as an ex- 
tended Slater-Pauling rule 

M tot = N V - 24n si (7) 

where nsi is the number of Si atoms in the unit cell. The 
Mn-rich structures show similar orbital-resolved DOS as 
in Fig. Mn atoms at Co sites mimic the DOS of the 
Co atoms in the minority-spin channel. Therefore, the 
original hybridization scheme proposed in Ref. [H is still 
valid and the SP rule also holds for the large unit cells of 
these compounds. The factor nsi takes the enlarged size 
of the unit cell into account. 

The CE for the total spin moment allows us to check 
millions of unknown structures whether they fulfill the 
extended SP rule. The result is plotted in Fig. [5] where 
the formation energies of structures that satisfy eq. ([7]) 
within a tolerance of 0.05 /is /atom (given by the cross- 
validation score) are highlighted. A large region of po- 
tential half-metallic compositions is identified. 

E. Monte Carlo simulations 



lation implied by the Slater-Pauling rule between the 
number of valence electrons and the total spin moment. 
These data are summarized in Tabic UTTl The spin gap de- 
creases with increasing Mn ratio until it is finally closed 
at a Mn concentration of about 60%. Thus, there is a 
wide Mn concentration range from 25% to 60% where 
the half-metallicity is preserved, and the total spin mo- 
ment per unit cell remains to be an integer multiple of 



From our defect calculations it has become clear that 
even ideally stoichiometric Co2MnSi shows a spin gap 
only if it is highly ordered. Hence the question arises 
how the properties of Co2MnSi are affected by thermal 
disorder at finite temperature. The ECE offers an excel- 
lent starting point to address this question by means of 
Monte Carlo simulations. The thermodynamic quantity 
governing the behavior of the material at finite temper- 
ature is the free energy, F = E — TS. While configu- 
rational entropy is treated implicitly in the Monte Carlo 
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simulation, the contributions from vibrational or mag- 
netic excitations should be included in the free energy 
functional employed in performing the MC simulations. 

As a first example, we restrict ourselves to CoaMnSi 
with ideal stoichiometry. In this case, thermal disorder 
will manifest itself by atomic swaps where we restrict 
ourselves to Co-Mn swaps that are detrimental to the 
half-metallic property. For a Monte Carlo simulation of 
stoichiometric Co2MnSi at finite temperature, the crucial 
input data are the free energy differences associated with 
defect formation. Analyzing the role of vibrational and 
magnetic excitations, we will conclude that none of them 
contributes significantly to free energy differences. Con- 
cerning the vibrational contributions, we argue that both 
the mass and the bond strength of Mn and Co atoms 
are very similar. Therefore, we expect only very small 
changes in the phonon density of states between ideal 
and defective samples. The phonon contribution to free 
energy differences is therefore neglected. As far as mag- 
netic excitations are concerned, we note that most defects 
couple ferromagnetically to the (ferromagnetic) host ma- 
terial. Far below the Curie temperature of Tq — 985 K, 
the magnetic excitations will consist of small oscillations 
around the direction of the average magnetization, both 
for the ideal and the defected materials. Because of this 
similarity, the magnetic contributions to the free energy 
for both situations tend to cancel. The only defect whose 
magnetic moment aligns antiparallel to the host material 
is the Mn atom on a Co site. For this case, we calcu- 
lated the energy to invert the Mn magnetic moment to 
be E mag = 1.1 eV (additional to the structural formation 
energy). Hence the Mn magnetic moment cannot flip 
freely, but, also in this case, will perform only small os- 
cillations around the ground state. Therefore, we neglect 
the magnetic contributions to the free energy differences 
as well. 

We perform Monte Carlo (MC) simulations of stoichio- 
metric Co2MnSi to sample the configurational space in 
thermodynamic equilibrium. Starting from a random oc- 
cupation of the lattice, the Mn atoms swap their sites 
with nearest- or next-nearest Co neighbors. The energies 
of the configurations are computed with the ECE Hamil- 
tonian, and the new configuration is accepted according 
to the Metropolis criterion. The MC results presented 
here have been obtained in a simulation cell containing 
6144 sites for Co and Mn atoms. The simulations show 
that the concentration of Co-Mn swaps is very low, below 
1 %o at temperatures of up to 1000 K (Fig. [9] , diamond 
symbols). Tests with an increased cell size of 12,000 sites 
gave the same results. We note that the interaction en- 
ergy between Mn atoms at Co sites is positive, but small 
(1.7 meV/atom from DFT-GGA, 2.9 meV/atom from 
ECE). At the low defect concentrations observed, the 
Co-Mn swaps can be treated as statistically independent 
and their concentration is given by 
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FIG. 9: (color online) The temperature dependence of the 
fraction of Co atoms at Mn sites S for different composi- 
tions in the simulation cell (blue circles: 4144 Co and 2000 
Mn atoms, red squares: 4120 Co and 2024 Mn atoms, black 
diamonds: 4096 Co and 2048 Mn atoms (stoichiometric com- 
position), green triangles: 4072 Co and 2072 Mn atoms, violet 
crosses: 4048 Co and 2096 Mn atoms). At temperatures be- 
low 1000 K the stoichiometric sample shows only l%o of Co 
defects and is almost perfectly ordered. Even above 1000 K, 
the Co antisite formation may be suppressed by excess Mn 
atoms. 



The comparison of the MC results with the analytical 
solution shows excellent agreement, confirming that it is 
indeed justified to neglect defect interactions in this case. 

Next, MC simulations are performed for samples that 
are slightly non-stoichiometric. We define the relative 
fraction S of the number of Co atoms nc at Mn sites 

Nmh, 



S = 



N M n 



(9) 



1 

CCo-Mn = — ^ exp 



E 



/ 



2k B T 



(8) 



Again, Fig. [5] shows that the system tends to an almost 
perfectly ordered state for temperatures below 1000 K. 
For higher temperatures, some additional Co antisites are 
created due to thermal disorder. For Mn-rich composi- 
tions, we note that the increase of thermally created Co 
defects is somewhat smaller than for the stoichiometric 
sample, i.e., the presence of excess Mn allows the sup- 
pression of Co defect formation. Moreover, we observe 
that the positive interaction energy between Mn atoms 
at Co sites leads to clustering of Mn antisites. 

We also perform MC simulations for a sample whose 
composition matches the newly predicted ground state 
Co2Mn4Si2 in order to assess the thermal stability of 
this ordered structure. First, we note that the en- 
ergy difference between the formation energy of ground 
state Co2Mn4Si2 and a linear interpolation between the 
energies of ideal Co2MnSi and Mi^Si is very small 
(5.6 meV/atom). The MC simulation shows that this 
leads to a phase separation of Co2Mn4Si2 into Co2MnSi 
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FIG. 10: Energy per atom vs. temperature from a MC sim- 
ulation with 4000 Co and 8000 Mn atoms in the cell cor- 
responding to ground state Co2Mn4Si2 at T = K. The 
dashed line shows the formation energy difference between the 
ground state Co2Mn4Si2 and a linear interpolation between 
ideal Co2MnSi and MnsSi. The intersection point coincides 
with the maximum slope of cy = 8E/dT\v that indicates the 
latent heat of the phase transition from the ordered ground 
state to phase separation at T « 320 K. 

and Mi^Si in thermodynamic equilibrium for tempera- 
tures T > 320 K (Fig. HDD. Thus, the potentially half- 
metallic ordered structures that we predict in a wide com- 
position range may not be practically relevant. Rather, 
the samples in this composition range could be described 
by precipitates of MnaSi (which is not a half metal) in a 
matrix of half-metallic Co2MnSi. The electrical proper- 
ties will then depend on the size and the concentration 
of these precipitates and on their percolation length. 

IV. CONCLUSIONS 



Both analytical calculations and Monte Carlo simula- 
tions for Co2MnSi close to the stoichiometric composition 
show that its thermodynamic ground state is almost per- 
fectly ordered up to temperatures of about 1000 K, with 
a Co-Mn swap concentration below l%o in thermody- 
namic equilibrium. The magnetization and especially the 
spin polarization at finite temperatures have been inves- 
tigated with an extended Heisenberg model (mapped to 
exchange interactions obtained from DFT calculations) 
and Monte Carlo simulations. It has been argued that 
the spin polarization of ordered Co2MnSi is nearly 100% 
for T < 0.27 T c (265 K)42 Therefore, it should be pos- 
sible to prepare samples (e.g. by careful annealing) that 
have a high structural order (avoiding the defect states) 
and indeed display the theoretically predicted spin gap 
at low temperatures. Adding more Mn to Co2MnSi can 
be used to keep the Co antisite concentration low even at 
higher temperatures. For non-stoichiometric, moderately 
Mn-rich Co2MnSi we have shown that the deficiency of 
electrons due to replacing Co atoms by Mn atoms is fully 
accommodated in the majority spin channel, and thus 
does not affect the half-metallic gap (extended Slater- 
Pauling rule). Using the cluster expansion for the total 
spin moment, we identify a large interval of Mn concen- 
trations ranging from 25% to 60% where potentially half- 
metallic compositions may occur. 

At low temperature, the cluster expansion for the total 
energy predicts four new ground state structures besides 
the well-known L2i structure. These new ordered struc- 
tures are slightly lower in energy than a random occupa- 
tion of Mn and Co sites. One of these new ground states, 
Co2Mn 4 Si2, even shows a spin gap. However, Monte 
Carlo simulations confirm the expectation that the order- 
ing energy associated with these structures is too small 
to render them stable at room temperature. Instead they 
decompose into regions of Co2MnSi and Mn 3 Si. 



In this work, cluster expansions on the basis of DFT- 
GGA calculations are reported for both the formation 
energy and the total spin moment of the pseudo-binary 
Heusler alloy Co2- x Mni+ x Si. These cluster expansions 
have been employed to study the thermodynamics of this 
alloy and its magnetic properties at various compositions. 
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